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We examine the threshold of the influence of the tunneling nanotubes (TNTs) on the cell-to-cell 
communication integrity. A deterministic model is introduced with the Michaelis-Menten dynamics 
and the intercellular exchange of substance. The influence of TNTs are considered as a functional 
perturbation of the main communication and treated as the matrix nearness problems. We analyze 
communication integrity in terms of the pseudospectra of the exchange, to find the distance to 
instability. The threshold of TNTs influence is computed for Newman-Gastner and Erdos-Renyi 
gap junction (GJ) networks. 
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Maintaining the functional integrity of cell-to-cell communication in multicellular systems is one of the prerequi¬ 
sites for achieving functionality in them m- In eukaryotic cells, intercellular communication is primarily mediated 
locally through gap junctions (GJs) and synapses; however, recent reports demonstrate the existence of a network 
of intercellular membrane nanotubes enabling long-distance communication EHH]. Also in prokaryotes they enable 
interspecies communication and share of antibiotic resistance liiin]. It has been shown that in both cases these tun¬ 
neling nanotubes (TNTs) can facilitate cell-to-cell communication and intercellular transfer of cytoplasmic molecules, 
organelles and viruses lainHn]. Existence of clusters of TNTs enables formation of complex cellular networks with 
both local and long-distance interactions based on membrane continuity between TNT-connected cells (Fig. [^. Em¬ 
pirical evidence indicate that they can have important role in many pathophysiological processes, like in activation of 
natural killer cells, regulation of osteoclastogenesis or in tumor formation and growth [IMS]. In prokaryotes they can 
play important part in transferring virulence from pathogenic to non-pathogenic bacteria |S] These findings indicate 
the need to systematically explore how perturbations in communication, induced by existence of clusters of TNTs, 
can influence integrity of intercellular communication. 

In this Letter we explore how the substance exchange through TNTs affects the functional stability of a multicellular 
system. We focus on two issues: (i) whether TNTs can either stabilize or destabilize intercellular communication 
governed by GJs? and (ii) how to determine the threshold of TNTs influence that shifts stable GJs communication to 
the unstable one? Therefore, we introduce a model of the substance exchange in a multicellular system, represented 
by ordinary differential equations, where cell-to-cell communication is mediated by both the GJs and TNTs while 
metabolic processes in cell follow Michaelis-Menten dynamics. In this model the GJs function governs the time 
evolution of the intercellular network while the TNTs function simulates the exchange mediated by the TNTs including 
a scaling parameter for that mediation. We consider the influence of TNTs as a functional perturbation of the main 
communication mediated by GJs. To determine the threshold for which the multicellular system remains stable, 
despite TNTs influence, we compute the pseudospectra of the exchange to find the distance to instability m- 

General model dynamies and pseudospectra. To investigate how TNTs affect the stability of the intercellular 
communication network, we model network dynamics as: 


x{t) = 'it{x{t)) := <^{x{t)) + ^e{x{t)), 


( 1 ) 


where x = [xi,X 2 , ■ ■ ■ tXuY". Here Xi{t) G [0,1] is the concentration of molecules and ions in the cell i £ N := 
[1, 2,..., n], d) : i?" —)■ ii" is a GJs function that governs time evolution of the intercellular network, while 0 : i?" —>■ i?" 
models the exchange mediated by TNTs and ^ > 0 is a scaling parameter for that mediation. TNTs are dynamic 
structures whose formation and decomposition is very sensitive to both intra- and extracellular factors, with lifetimes 
ranging from several minutes up to 4-5 hours iHiiiaiisj. Their short lifetime suggests that they can promote unstable, 
transient cell-to-cell communication, in contrast to more stable GJs HHj. Since many questions remain unanswered 
about how cargo is transported through TNTs [T9| we consider their effect on the model dynamics as an uncertainty 
described by 0. The system 0 is generally a nonlinear one whose stability is typically investigated at the equilibrium 
state X* G i?" as the local asymptotic stability, where x* is a state vector such that 4'(a:*) = 0. An equilibrium state 
X* is locally asymptotically stable if there exists e > 0 such that for every a:(0) that is in e neighbourhood of the 
equilibrium x* (i.e., jja;* — x(0)|| < /r), it holds that lim \\x{t) — x*|| = 0. This stability property is characterised by 


the spectra of the Jacobian matrix A = 



of (1) in the state x = x*, which can be written as A = A -|- ^A, 



FIG. 1. Intercellular TNTs between neighboring cells. A field of cells with a cluster of smaller TNTs (highlighted by a dashed 
circle) and a more pronounced larger tube (indicated by an arrow). Reprinted with permission from [Sj. 
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where A = corresponds to the exchange through GJs that will be named as the observable part. The 

term A = is determined by the transport through TNTs that is generally unknown. Thus, we call it 

the unobservable part. Furthermore, for ^ > 0, to reflect the scale of TNT mediation we assume that uncertainty 
parameters are unit scaled in the chosen matrix norm, i.e. ||A|| = 1. 

To investigate how TNTs can influence stability of communication we will analyze sensitivity of the spectrum of 
the observable matrix A upon perturbation ^A, using the concept of matrix pseudospectra |20j . Given a matrix 
A e C"’” and e > 0, the e -pseudospectrum of a matrix A, denoted by A£(^), is composed of all eigenvalues of 
matrices which are ”£—close” to A: A S Ae(A) if and only if there exists x G C" \ 0 and A G C”’” such that |j A|| < e 
and (A + A)x = Xx, i.e. Ai;{A) = U||A||<e^(^ + ^)- Consequently, we use e-pseudospectrum to establish spectral 
properties that are robust under matrix perturbations bounded in a given norm || • || by the parameter e > 0. The 
measure of robustness of the stability of the observable matrix A is defined as the largest £ > 0 such that Ae(A) C C“, 
and is known as the distance to instability m- Since the computation of such value £ > 0 requires solving a nonconvex 
optimization problem, numerical algorithms have to be used pnil24) . Computing the spectra A (A) of the observable 
part of the Jacobian, corresponding to GJs interactions and intracellular metabolic processes, we can determine the 
expected stability A(A) C C“ or instability A(A) ^ of the substance fluxes. However, the unobservable part 
of cellular communication can change this spectral property and lead to the different dynamics of the network since 
A(A) C A^(A). More precisely, assuming that A(A) C C“, then A(A) C Aj(A) C C“ if and only if C < where zzZ 

denotes the distance to instability of A. Therefore, if the TNTs influence is dominated by the distance to instability of 
the observable matrix A, we can safely conclude that the mediation of TNTs does not change the stability of Q, while 
for the values ^ > ly^we cannot do that. Otherwise, assuming that A(A) % C“, then A(A) % C“ if and only if ^ 
where izi denotes the distance to stability of A. Thus, in this case, distance to stability of the observable matrix 

A is the critical value of the TNTs influence for which the system can change the unstable behaviour to the stable 
one. Therefore, the integrity of GJ communication (either stable or unstable one) under the influence of TNTs does 
not depend only on GJ network’s resilience^ i.e. real part of the least negative eigenvalue of the observable Jacobian 
matrix. In fact, it remarkably depends on the GJ network structure which can produce non-normal dynamics |2()j . 

Deterministic model. We introduce a deterministic model for substance exchange in a multicellular system which is 
mediated by GJs and TNTs. Both, GJs and TNTs allow various molecules and ions to pass freely between cells through 
the channels by the diffusion like process. However, diffusion in cells and between them, known as anomalous diffusion, 
can differ from classical one due to spatial inhomogeneity p51127) . In situation like this, it is suitable to consider the 
kinetics of substance exchange between cells in terms of an exchange coefficient gij with dimensional unit of inverse 
time |28] . In the simplest case, communication from the cell i to cell j is proportional to the concentrations between 
the cells. Therefore, we can define the substance exchange between cell j to cell i as gij{xj{t) — Xi{t))+^5ijXj{t), where 
^ > 0 is a small value that determines the strength of influence of TNTs on communication modeled by uncertainty 
parameter 5ij . Since exchanged substances play a role in the metabolic processes inside the cell and are released into 
environment, we introduce parameters at > 0 that describe the rate by which the substance is metabolized by the 
cell i G N in time t. Since most of metabolic processes follow Michaelis-Menten dynamics, we introduce /3i > 0 as the 
half-time saturation coefficient for the cell i. Therefore, we express the intercellular communication as: 


Xi{t) 


OliXiit) 

+ Xi{t) 


+ XI 3i3 + C X (* € A^), 




jeN 


( 2 ) 


Finally, we restrict to the case of zero equilibrium x* = 0, when the initial state a;(0) A 0 reflects starting distribution 
of the substance in the network. Gonsequently, A = [J^-] with ||A|j = 1, while A = [§|^(0)] = [a^-] where for G N 
dij = — (^ + J2j^i9ij) if * = J 01' = 9ij if * A j- We use such a simple deterministic model because: (i) it is 

general enough to embrace the main aspects of cell-to-cell substance exchange via GJs and TNTs and (ii) it, without 
the loss of generality, clearly illustrates the use of pseudospectra in estimating the uncertainty in the model dynamics 
determined by the substance exchange through TNTs. 

Depending of the system’s quantity or property which is examined, we chose the norm in which distances are 
measured. Here we discus two cases. First, that is the maximal substance concentration in the network with norm 
set to infinity norm. In that setting, for communication mediated exclusively by GJs, i.e. ^ = 0, using pseudospectral 
localization ^ , we can conclude that zero equilibrium point of (H) is exponentially stable one, and that for all ^ A 0, 
the flow of x{t) satisfies ||a;(t)||oo A e“^+'’||a;(0)||oo for t > 0, where M = minigA? ^ > 0. In other words, for 
f < M, the maximal substance concentration in individual cells exhibits an exponential decay without the initial 




4 



0.025 

0.02 

0.015 

0.01 

0.005 

0 

- 0.005 

- 0.01 

- 0.015 

- 0.02 

- 0.025 



time [sec] 


FIG. 2. Distance to instability computed for 100-cell GJ network constructed as the Newman-Gastner spatial network, (a) 
Graphical image of that network; (b) pseudospectral portrait of its Jacobian matrix A (asterisks mark eigenvalues of (^ = 0), 
while the shadowed area represents all the possible locations of eigenvalues for the full network Jacobian A, when GJs and 
TNTs are included = ^crit where ^crit = is a threshold value)); (c) transient growth of substance concentration within 
cells from the initial state measured in the Euclidian norm ||e^^|| 2 , due to non-normality of the GJs Jacobian matrix |34] for 
the following cases: ^ = 0 (solid black), ^ = ^crit (solid gray), ^ = O.S^cHt (solid dashed-dotted) and ^ = 1.2^crit (solid dashed). 
Note that peak of concentration and duration of concentration decay depends on the scenario used as seen in Figs 2c-4c. The 
computed values are ^crit = 2.38 • 10■^ and M = 4.47 • 10“®. 


growth. This implies that the system maintains its communication integrity. Therefore, the obtained constant M > 0 
is a lower bound of the distance to instability of >1 in || • ||oo. Second, to check the behavior of the network in terms of 
total square deviation from the equilibrium state, the norm has to be set to the Euclidian one. Therefore, to compute 
the threshold of TNTs influence, in that setting, we apply the algorithm developed by Kostic et al. m, while the 
pseudospectral portraits are generated using the EigTool [3T]. 

Results and discussion. We consider a few realistic scenarios based on the simple deterministic model of 100-cell 
network. In all test cases we compute lower bound of threshold of TNTs influence in infinity norm (M), the exact 
value in Euclidean norm and construct pseudospectral portrait with transient plot. Eor GJ communication we 
use: (i) spatially distributed Newman-Gastner network with weight parameter set to 0.001 (35] and (ii) Erdos-Renyi 
modular network |33j with ten clusters connected wit 0.03 overall probability of the attachment and proportion of 
links within modules set to 90%, graphically depicted in Figs. and For values of physiological parameters 

in simulations, we assume that all cells in the population are of the same type and therefore have the same time scales 
for metabolizing the substance. Therefore, for each i = 1,..., 100, we take Pi randomly on an uniformly distributed 
interval [0.9,1.1]. On the other hand, the saturation constants (having the same order of value as the exchange 
coefficients), that are of the same order as the exchange coefficients, can differ more significantly, and are chosen 
randomly from [0,0.01] with the uniform distribution. Finally, following ^7\, for exchange coefficients gij we use 
random values from the interval [0, 0.05] with the uniform distribution. 

We compute the critical pseudospectra of the GJs Jacobian matrix A of the network (Figs. [^,[^ andj^)). Here, 
the term critical pseudospectra stands for the fact that e = which is the computed threshold of TNTs influence. 
The shadowed area indicates how the system is sensitive to changes in cell communication determined by TNTs. To 
demonstrate how formation of TNTs can affect network dynamics, we compute first order approximated behavior of 
([|), measured in Euclidean norm, i.e. Ile*^ll 2 (Figs. and|^), for the following cases. Case 1, when cell-to-cell 

communication takes place only through GJs (^ = 0). Then the system, after passing through short transient interval, 
is reaching the stability either faster (Figs. and[^) or slower (Figs. |^). The corresponding curves are depicted 
by the solid black lines. Case 2 (critical case), when ^ = ^crit and A is the worst arrangement of TNTs that move 
eigenvalues of A to the imaginary line {marginal instability). Such A is constructed via suitable singular vectors [30j . 
From Figs. and 1^ (solid gray line) it is seen that the substance exchange in the system is in the state of an 

oscillating mode, waiting to start towards either stability or instability. Case 3 (case of stability) when ^ < ^crit 
(here ^ = O.S^crit)- In this case the system maintains the communication integrity (gray dashed-dotted line in Figs. 

and|^). Case 4 (case of instability), when ^ > ^crit (here ^ = 1.2^crit)- Gorrespondingly, the system is 
communicationally disintegrated (gray dashed line in the Figs. and|^). 
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FIG. 3. The same as in Fig. but for Erdos-Renyi modular network. The computed values are ^crit = 3.15-10 and 
M = 6.30 • 10“®. 
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FIG. 4. The same as in Fig. [^but for Erdds-Renyi random network simulating a pathological state of the system (colored 
circles show nodes with altered capacity to receive signals). The computed values are ^crit = 6.55 • 10“^, and M = 6.30 • 10“^. 


The situation of disrupted cell-to-cell communication can cause numerous diseases. For example, inborn cardiac 
diseases are among the most frequent congenital anomalies and are caused by mutations in genes that form GJs |35) . 
Therefore, it is crucial to investigate stability of intercellular communication and determine possible thresholds for 
disruption of cell-to-cell communication integrity. To investigate a possible influence of TNTs in the case of disrupted 
communication we simulate pathological situation by modifying the Erdds-Renyi modular network as follows. In this 
network arrangement only one module (colored circle in Fig. exhibits a cascade degradation of their capacity to 
receive the substance under exchange, while their capacity to send it in the fixed network flux direction (corresponding 
to the node enumeration) remains the same. More precisely, in the original realization of Erdds-Renyi modular network 
we take O.lgtj, instead of gij for i = 1... 100 and j = 1... 10. In the example we create, when only one module is 
disrupted, pseudospectral portrait (Figj^) shows that sensitivity of the whole network to communication changes is 
significantly increased, compared to non-pathological state (Figj^). Also, critical level of oscillations in Euclidian 
norm deviate more from the equilibrium state indicating that formation of TNTs can disrupt the system more easily. 

To summarise, the main novelty of this Letter lies in the fact that the uncertainty of TNTs influence to the 
overall cellular communication can be treated as the matrix nearness problems, i.e. either as distance to instability 
(treated in this Letter) or distance to stability (reverse problem that will be treated in the follow-up study). We have 
presented how this concept can provide meaningful insights using the simple deterministic model of cellular networks 
with asymptotically stable GJs cell-to-cell communication, where TNTs can destabilize the system. The problem 
is analysed in terms of maximal individual deviation (|| • ||oo) and total square deviation (Euclidian norm) of the 
cells substance concentration. In the simulations, the threshold of such TNTs influence is computed using recently 
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developed pseudospectra methods for two standard structures of cellular networks (spatial Newman-Gastner and 
modular Erdos-Renyi) modelling healthy and pathological state of the system. The work presented here is a first step 
towards understanding of the influence of TNTs as uncertainty of the system using matrix analysis and computational 
methods. Therefore, many open questions remain, targeting, among other ones, the tie between matrix perturbations 
as stochastic processes, where the pseudospectra combined with the Bregman divergences [3S] can help to reliably 
estimate the mathematical expectation of the threshold of destabilising/stabilizing intercellular communication. 

This work is supported by the Ministry of Education and Science of the Republic of Serbia (Grants III 43007 and 
174019) and Provincial Secretariat for Science and Technology Development of Vojvodina (Grants 1136, 1850 and 
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